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PACS 47 . 35 . -i - Hydrodynamic waves 

PACS 68 . 03 . -g - Gas-liquid and vacuum-liquid interfaces 

Abstract. - We study theoretically the capillary-gravity waves created at the water-air inter- 
face by a small two-dimensional perturbation in the frequently encountered case where a depth- 
dependent current is present in the fluid. Assuming linear wave theory, we derive a general 
expression of the wave resistance experienced by the perturbation as a function of the current 
profile in the case of an inviscid fluid. We then illustrate the use of this expression in the case of 
constant vorticity. 



Water waves are both fascinating and of great practi- 
cal importance [l}|3]. For these reasons, they have at- 
tracted the attention of scientists for many centuries [I]. 
Water waves can for instance be generated by the wind 
at sea, by a moving boat on a calm lake, or simply by 
throwing a pebble into a pond. Their propagation at the 
surface of water is driven by a balance between the liq- 
uid inertia and its tendency, under the action of grav- 
ity or of surface tension (or a combination of both in the 
case of capillary-gravity waves), to return to a state of 
stable equilibrium [5]. Neglecting the viscosity of water, 
the dispersion relation of linear capillary-gravity waves 
relating the angular frequency u) to the wavenumber k 
is given by u 2 = (gk + jk 3 / p) tanh(kh), where 7 is the 
liquid-air surface tension, p the liquid density, g the ac- 
celeration due to gravity and h the depth of water [11. 
The above equation may also be written as a depen- 
dence of the phase velocity c = cu/k on the wavenumber: 
c(k) = (g/k + ^/k/p) 1 ^ 2 (tanh(fc/i)) 1 ^ 2 . The dispersive na- 
ture of capillary-gravity waves is responsible for the com- 
plicated wave pattern generated at the free surface of a 
still liquid by a moving disturbance such as a partially 
immersed object (e.g. a boat or an insect) or an external 
surface pressure source. The propagating waves generated 
by the moving disturbance continuously remove energy to 
infinity. Consequently, the disturbance will experience a 
drag, R, called the wave resistance j6j. In the case of boats 
and large ships, this drag is known to be a major source 
of resistance and important efforts have been devoted to 
the design of hulls minimizing it [7]. The case of objects 



small relative to the capillary length kT 1 — (7 /{pg)) 1 ^ 2 
has only recently been considered [8 ■ 14 and has attracted 
strong interest in the context of insect locomotion on wa- 
ter surfaces 15 16 . In the case of a pressure distribution 



of amplitude p localized along a line and traveling over 
the surface with speed V perpendicularly to its length, 
the wave resistance for deep water (h — > +00) is given by 



R = (p 2 /7)(l 
for V < c min 



(c min /vy 



for V > c„ 



and R = 



2,8 



Here c m i n = (^gj/p) 1 / 4 is the minimum 
of the wave velocity c(k) for deep-water capillary-gravity 
waves. In the limit V 3> c m j n , the wave resistance reduces 
to p 2 /"f. Note that as V approaches c m i n (from above), 
the wave resistance becomes unbounded. This diverging 
behavior of the wave resistance (in the case of a two- 
dimensional pressure distribution) is related to the fact 
that when V approaches c mm , the phase velocity ui(k)/k 
(equal to V, see [8j) and the group velocity dui(k)/dk tend 
towards the same value. It was shown in M that in the 
presence of viscosity, the wave resistance remains bounded 
as V approaches c m ; n . A similar regularization also exists 



when one takes into account non-linear effects 17 



In many cases of physical interest (like wind-generated 
flows 18j 19]), the waves propagate on shear currents 
rather than in still water (see, e.g., the seminal work of 
Miles (20] )■ The general problem of the interaction be- 
tween water waves and arbitrary steady current is of great 
physical significance 21 . It is, however, rather difficult 
and remains largely unsolved 22-26 . In this letter, we 



consider how the above predictions for the wave resistance 
are modified when the pressure distribution propagates on 
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Fig. 1: Schematic diagram: an external surface pressure dis- 
tribution moves with speed V in the x direction above the free 
surface of a liquid of depth h in which a velocity profile U(z) 
preexists. 



steady shear currents (taking also into account the effect 
of water finite depth, assuming h > \/3/2 

The present letter is organized as follows. We first for- 
mulate the problem by analyzing the fluid equation of mo- 
tion together with the boundary conditions. This allows 
us to get the free surface displacement as a function of 
the shear current and the pressure disturbance. We then 
establish a general expression for the wave resistance ex- 
perienced by the perturbation as a function of the current 
profile. Finally, we illustrate the use of this expression in 
the case of constant vorticity. 

Model - We consider the two-dimensional motion of a 
layer of fluid (assuming the fluid to be incompressible and 
inviscid). This implies that all physical quantities depend 
spatially on only one horizontal coordinate denoted by x, 
and on the vertical coordinate denoted by z. The flat 
bottom is given by z = — h while the free surface (in the 
absence of waves, see below) corresponds to z = 0. We 
assume the existence of a steady shear current below the 
free surface characterized by a velocity component U(z) 
in the horizontal x-direction, and a velocity component 
equal to zero in the vertical z-direction. In addition to 
this flow, capillary-gravity waves are generated by a pres- 
sure distribution (invariant along the ^-direction) moving 
with a constant speed V along the x-direction, as indi- 
cated in FigJT] Let z — C( x jt) denote the displacement 
of the free surface (in the presence of waves) . The veloc- 
ity component in the horizontal a;-direction is now given 
by U(z) + u(x,z,t), and the velocity component in the 
vertical z-direction by w(x,z, t). The additional velocity 
field (u, w) is assumed to be a first order correction to the 
undisturbed flow. Associated with the wave-induced mo- 
tion is a stream function ip(x,z,t) so that u = d z ip and 
w — —d x tp. Note that the use of the stream function does 
not impose an irrotational fluid motion. 

Having in mind that the phase velocity is imposed by 
the velocity V of the pressure disturbance [Tl[8l , we shall 



seek a stream function of the form 

dk 



i/)(x, z,t) 



2tt 



A(k)f k (z) e ik ^- vt l 



(1) 



We shall now determine the function f k (z) (characterizing 
the z-dependence of the stream function) using the equa- 
tion of motion and the boundary conditions. According 
to Euler's equation we have 

p(d t u+U(z)d x u + wU'{z)) = -d x p, (2) 
p(d t w + U(z)d x w) = -d z p~g. (3) 

Using ([!]), the above two equations can be rewritten as 
d k 

-A(fc) ikp{(V - U)f' k + U'f k }e ik ^~ v ^ = d xP , (4) 
dk 

^A(k) pk 2 {{V - U)f k } e^- y ^ = d z p + pg. (5) 

Eliminating pressure between equations Q and (|5| yields 

(V - U(z)) [ft{z) - k 2 f k (z)} + U"(z)f k (z) = 0. (6) 

Equation ([6]), which relates the function fk(z) to the cur- 
rent profile U(z), is known as the inviscid Orr-Sommerfeld 
or Rayleigh equation 22 26;. It has to be supplemented 



with boundary conditions. At the bottom, the fluid veloc- 
ity vanishes and so f k (z = —h) = 0. Let Uq = U(z = 0) 
and Uq = d z U(z = 0). The dynamic free surface boundary 
condition [2] can be written as 



dk 
2^ 



A{k) p{~(U - V)f' k (0) + U o f k (0)} e ik ^- vt ^ 

= p + pgC, (7) 



where, according to Laplace's formula, the pressure p 
equals p ext (x,t) — jd x C pi- Note that /fc(0) and f' k {0) 
depend on the full profile U(z) with z £ [—^,0] (see 
(g). Let Pext (x,t) = (27T)- 1 / dk P ext (k) e ^-vt) and 
C(x,t) = (27T)- 1 / dk((k)e lk( - x ~ vt l Using then the kine- 
matic free surface boundary condition w = d t Q + Uq8 x Q 
we obtain A(fc) = -({/ - U)C(fc)// fe (0) and 



p\ (U -V) 



a /fc(0) 
/fc(0) 



%{U -V) 

= Pext{k). 



C(fc) 



(8) 



Equation ([8]) is of physical importance since for a given 
current profile it relates the Fourier component £(&) of the 
surface displacement to the Fourier component P ex t(k) of 
the pressure disturbance. 

Let us emphasize that despite the fact that we are work- 
ing within the frame of a linear wave theory, a linear com- 
bination of solutions corresponding to different current 
profiles will generally not be a solution of the problem 



at hand. Sec 22 26 for further details on this issue. 

Wave resistance - We can now investigate the wave 
resistance R experienced by the disturbance. According 



p-2 



Capillary-Gravity Waves on Depth-Dependent Currents 



to Havelock [6], we may imagine a rigid cover fitting the 
surface everywhere. The pressure p ex t{x) is applied to the 
liquid surface by means of this cover; hence the wave resis- 
tance is simply the total resolved pressure in the x direc- 
tion. This leads to R — — f dx p(x) d x ((x) |6j. According 
to ^ the wave resistance can then be written as 



R 



dk 



Ak\P ext (k){< 



^iuo-vym-u^-v)- 



P 



(9) 



which is the central result of the present letter. It allows 
one to calculate the wave resistance experienced by the 
moving disturbance for any current profile U{z). Note that 
if one replaces V by ui/k and sets the denominator of ([9| to 
zero, one gets a general formula for the dispersion relation 
of capillary-gravity waves on depth-dependent current. 




Fig. 2: (Color online) Plot of the wave resistance R* for 
V — in units of p 2 /2j, as a function of the reduced veloc- 
ity C/o/cmin- The red curve (dotted) corresponds to a uniform 
current U(z) = Uo with infinite depth, the black one (dashed) 
corresponds also to a uniform current U(z) = Uo with finite 
depth h — 5k _1 , while the blue one (solid line) corresponds to 
a linear current U(z) — Uo(h + z) /h with the same finite depth 
h = 5k _1 . The insert shows the corresponding dispersion rela- 
tions. 



The integral in equation ([9| cannot be evaluated unam- 
biguously because the poles of the integrand are on the 
domain of integration. This ambiguity can be removed by 
imposing the radiation condition that there be no wave 
coming in from infinity. There are several mathematical 
procedures equivalent to this radiation condition. One 
way is to consider that the amplitude of the disturbance 
has increased slowly to its present value in the interval 
— oo < t < : pl xt (x,t) — e et p ext (x,t) where e is a small 
positive number that will ultimately be allowed to tend to 
zero. The poles of the integrand have been shifted over 
and above the real axis. Since the poles are now out of 



the domain of integration, the integral can be evaluated 
numerically unambiguously. 

For a non linear current profile, i.e. U"(z) ^= 0, 
one cannot in general find explicit solutions to the Orr- 
Sommerfeld equation In that case the wave resis- 

tance has to be evaluated numerically. For further de- 
tails on approximation methods for the Orr-Sommerfeld 
equation, readers may refer to [22j|26]. In what follows 
we will thus assume the shear current U (z) to be linear: 
U(z) = Uo(h + z)/h, corresponding to a constant vor- 
ticity zu = —Uo/h [27]. This might be of relevance for 
tidal waves 23, 28] and constitutes a first step towards 
more general cases. In that case U"(z) = and the Orr- 
Sommerfeld equation (JoJ) admits an exact solution 



f k (z) = (V-U ) 



sinh(fc(z + h)) 
smh(kh) 



(10) 



Inserting (10 1 into ([9]), the wave resistance is then given 
by 



R = lira 



dk -ik\P &xt {k)f 



>o J 2np 

where the denominator (...) corresponds to 



(11) 



(U - Vf 



2i(U - V)e 
k 



U' [Uo - V - - 



k coth(kh) 

?+ — 
P 



(12) 



For a well localized pressure distribution of the form 
pix) = p8(x), one has P ex t(k) — p. The corresponding 
wave resistance is shown in Fig. 2 as a function of the re- 
duced velocity Uo/c m - ln assuming V — 0. To gain physical 
insight in the problem we also display in Fig. 2 the wave 
resistance for a uniform current U(z) = Uo, both for finite 
and infinite depth. We first observe, in the case of a uni- 
form current with finite depth (black dashed curve), the 
existence of a shoulder at Uo = y/gh- This is related to 
the fact that when k — > 0, the phase velocity uj/k tends to 
a finite value \fgh (see insert in Fig. 2) [2j. Hence, when 
Uo becomes larger than y/gh, one of the two poles in the 
integrand of equation ( [IT] ) disappears. At first sight it is 
therefore surprising that no such shoulder is present in the 
case of linear current U(z) = Uo(h + z)/h (blue solid line). 
In fact, one can show that in the case of linear current and 



increases 29 . 



V = the integrand of ( 11 1 conserves its two poles as U 



We may also emphasize the fact that for both finite 
depth cases with uniform and linear current the minimum 
phase velocity is shifted from its value c m i n = (^gj/p) 1 ^ 4 
defined in the introduction of this letter for infinite depth. 
Therefore, in Fig. 2, the solid blue curve as well as the 
dashed black curve do not diverge exactly at abscissa 
Uo/c m - m = 1 but for quantities slightly superior to it. Nev- 
ertheless, for h larger than say 3k _1 , the minimum phase 
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velocity for the uniform current and finite depth case be- 
comes very close to c min . 
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Fig. 3: (Color online) Plot of the free surface profile C,(x,t) in 
units of p/27, as a function space x for Uo/cmin = 1.2 and 
h = 5k" 1 . The red curve (dotted) corresponds to a uniform 
current with infinite depth, the black one (dashed) corresponds 
also to a uniform current with finite depth h, while the blue one 
(solid line) corresponds to a linear current U(z) = Uo(h + z)/h 
with finite depth h. 

Fore completeness, we have displayed in Figj3] the ele- 
vation of the free surface in the case of a linear current 
as obtained from Q (assuming a well localized pressure 
distribution of the form p(x) — pS(x)) and compared it 
with the surface elevation for uniform currents (finite and 
infinite depth). 

Conclusion - We have shown that non uniform currents 
have important effects on capillary-gravity waves gener- 
ated by a two-dimensional perturbation. Both the waves 
properties and the corresponding wave resistance are sig- 
nificantly modified when a current exists in the fluid, ft 
would be of great interest to extend the results of the 
present letter to three-dimensional perturbation. A deeper 
understanding of the physical response of the wave system 
near c m j n will also require the introduction of nonlinear ef- 
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